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Spin polarization and g-factor of a dilute GaAs two-dimensional electron system 
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(Dated: February 1, 2008) 

The effective g-factor (<?*) of a dilute interacting two-dimensional electron system is expected 
to increase with respect to its bare value as the density is lowered, and to eventually diverge as 
the system makes a transition to a ferromagnetic state. We report here measurements of g* in 
dilute (density 0.8 to 6.5 xl0 10 cm~ 2 ), high-mobility GaAs two-dimensional electrons from their 
spin polarization in a parallel magnetic field. The data reveal a surprising trend. While g* is indeed 
significantly enhanced with respect to the band g-factor of GaAs, the enhancement factor decreases 
from about 6 to 3 as the density is reduced. 

PACS numbers: 73.50.-h, 71.70.Ej, 73.43.Qt 
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The ground state of a dilute, interacting electron sys- 
tem has been of interest for decades. It has long been ex- 
pected that, because of interaction, such a system makes 
a transition to a ferromagnetic state as the density is re- 
duced below a certain threshold jl], |j. For even lower 
densities, the system should eventually become an elec- 
tron solid (Wigner crystal) |J . A relevant parameter as- 
sociated with this evolution is the effective g-factor (g*) 
of the system. In the limit of high density, when the 
parameter r s , the average interparticle separation mea- 
sured in units of effective Bohr radius, approaches zero, 
g* should have the " bare" value determined by the energy 
band structure of the host material. With decreasing 
density, g* is expected to increase monotonically and di- 
verge at the density below which the electron system en- 
ters its ferromagnetic state. For an ideal two-dimensional 
electron system (2DES), quantum Monte Carlo calcula- 
tions pL M indeed confirm the above trend. 

An excellent candidate for testing these predictions is 
the GaAs 2DES, as it possesses very low disorder com- 
bined with a simple band structure. Here we report mea- 
surements of the spin polarization of a very high quality, 
dilute 2DES in a modulation doped GaAs/AlGaAs het- 
erostructure as a function of an in-plane magnetic field. 
Via transport measurements, we find the magnetic field 
above which the 2DES becomes fully spin polarized, and 
from this field we determine g* . The results reveal a re- 
markable trend: as the density is lowered from 6.5 to 0.8 
xl0 10 cm~ 2 , corresponding to an increase in r s from 2.1 
to 6.3, the measured g* decreases from 2.7 to 1.3. This 
implies a substantial overall enhancement of g* with re- 
spect to the band g-factor of GaAs (\gb\ = 0.44 in GaAs 
||). The decrease of g* with r s , however, is unexpected. 

We studied a Si-modulation doped GaAs heterostruc- 
ture grown on a (100) GaAs substrate. We used a square 
sample in a Van der Pauw geometry, with a backgate to 
control the density. We made measurements in a dilu- 
tion refrigerator at a temperature (T) of ~ 25mK and 
magnetic fields (B) up to 18T, and in pumped 3 He at 
~ 0.3K and fields up to 33T. The sample was mounted 
on a single-axis tilting stage that can be rotated, using 
a computer controlled stepper motor, in order to change 
the angle (9) between the sample plane and the magnetic 



field. The measurements were done using low-frequency 
lock-in techniques. At zero gate bias, the sample has a 
density n = 1.4 x 10 10 cm~ 2 and a mobility of 55 m 2 /Vs. 
An external magnetic field applied parallel to the 2DES 
causes a Zeeman splitting of the energy bands. This 
splitting induces a difference in population of the spin- 
up and spin-down subbands, which leads to a net spin 
polarization of the system. If the splitting exceeds the 
Fermi energy of the system, all spins are aligned and the 
2DES is fully spin polarized. Assuming a simple model, 
in which g* is independent of the applied magnetic field, 
we can write the splitting between the spin-up and spin- 
down subbands as Ez = |<7*|Mb-B , where jib is the Bohr 
magneton. In this model the 2DES becomes fully spin 
polarized at a field Bp, given by 
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where m* is the effective mass, n is the total density of 
the 2DES and h is Planck's constant. We emphasize that 
we measure the effective g-factor defined, as in Eq. (1), 
by the field at which full polarization is achieved. 

In our experiments we measure, via the analysis of 
Shubnikov-de Haas (SdH) oscillations, the Fermi con- 
tours of the two spin subbands. We apply a constant 
magnetic field parallel to the 2D plane and slowly rotate 
the sample around 9 = 0° to induce a small perpendicu- 
lar field component, B±. We record the sample resistance 
during the rotation, and Fourier analyze its SdH oscilla- 
tions with B± to obtain the populations of the two spin 
subbands M . Our experiments allow a determination of 
the field, Bp, above which the minority spin-subband 
depopulates and the 2DES becomes fully spin polarized. 
In the range where the Fourier transforms are done the 
parallel component of the field, Sii , is equal to the total 
field B to better than 2%. 

In Fig. 1 we show plots of the sample resistance R vs. 
B±, taken at a density n = 2.05 x 10 10 cm~ 2 , as deter- 
mined from the positions of quantum Hall states. The 
top trace was taken in a purely perpendicular field. The 
Fourier transform (FT), shown on the right, exhibits two 
peaks, one at 0.85T and another at approximately half 
this value, 0.42T. The 0.85T frequency, when multiplied 
by (e/h), gives 2.05 x 10 10 cm -2 , i.e., the total density 
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FIG. 1: Resistance vs. perpendicular field at the indicated 
parallel fields and the Fourier transforms of the SdH oscilla- 
tions. The total density of the 2DES is 2.05 x 10 10 cm -2 . 



of the 2DES. The 0.42T peak stems from the spin unre- 
solved SdH oscillations. The rest of the traces shown in 
Fig. 1 were taken by rotating the sample at the indicated 
B applied almost parallel to the 2DES. With increasing 
Bn, we observe a splitting of the lower FT peak (0.42T) 
into two peaks. The positions of these two peaks, mul- 
tiplied by (e/h), give the two spin subband populations. 
Note that the two populations add up to the total den- 
sity of the sample. As B» is increased, the majority spin 
subband peak merges with the total density peak (0.85T) 
and the minority spin subband peak moves to very low 
frequencies and is no longer resolved ||. 

In Fig. 2 we summarize the positions of the FT peaks 
corresponding to the majority and minority spin sub- 
bands as a function of B for the case examined in Fig. 

1. Above a certain field Bp the majority spin subband 
population saturates at a value which corresponds to the 
total density of the 2DES. Therefore, Bp marks the onset 
of full spin polarization. Within the experimental error, 
the evolution of the FT peak positions as a function of 
field is linear. This implies that g* is roughly indepen- 
dent of the applied parallel field. 

We also measured the in-plane magnetoresistance 
(MR), by fixing 9 at 0° and recording the resistance as 
a function of the applied magnetic field. The MR trace, 
taken for n = 2.05 x 10 10 cm -2 , is also shown in Fig. 

2. This trace exhibits a clear break in the functional 
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Summary of FT peak positions as a function of ap- 
plied magnetic field, along with the in-plane magnetoresis- 
tance trace. The dashed line marks the field Bp at which the 
2DES becomes fully spin polarized. 



form of the MR as it changes from an ~ e B dependence 
at low field to a simple exponential, ~ e B dependence 
at higher fields. The data in Fig. 2 demonstrate that 
the onset of the simple exponential behavior of the in- 
plane MR coincides with the field Bp above which the 
spins are fully polarized. Remarkably, the same func- 
tional behavior of the in-plane MR is seen in GaAs 2D 
hole systems M, and it has been shown J7J that the on- 
set of the simple exponential regime corresponds to the 
full spin polarization of the 2D system. Several studies 
on 2D electrons in Si-MOSFETs have also pointed out 
a correlation between the in-plane MR and the full spin 
polarization [|l0|, [0], g |D|. 

The strong correlation between the full spin polariza- 
tion and the onset of exponential MR provides another 
method of finding the field Bp. While this method does 
not allow a direct measurement of the spin subband pop- 
ulations, it can be useful at lower densities where the 
SdH method is no longer practical because of a decrease 
in the number of resistance oscillations. We summarize 
in Fig. 3 our MR data taken at different densities. At all 



densities, the MR exhibits an 



dependence at high 



fields and the onset of this dependence clearly depends 
on the density. We emphasize that in several cases where 
we have made both MR and SdH oscillations measure- 
ments in constant field, the values of Bp obtained from 
the two methods coincide. 

Using the measured values of Bp and the relation 
(I) we determine g*. For Bp we use the saturation 
field of the majority spin subband FT peak and the on- 
set of the exponential regime of the in-plane MR. An- 
other parameter needed is the effective mass, to*, which 
also can be different from the band effective mass due 
to electron-electron interaction |L4j. We independently 
measured to* from the T dependence of SdH oscilla- 
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FIG. 3: In-plane magneto-resistance at different densities, 
from top to bottom: 0.8, 1.08, 1.4, 1.68, 1.9, 2.18, 2.47, 
2.89xl0 10 cm -2 . Dashed lines mark the field above which 
the magnetoresistance exhibits a simple exponential depen- 
dence. Top two traces were taken at a temperature of 0.3K, 
the rest at 25mK. Top trace was divided by a factor of 3. 



tions in a purely perpendicular magnetic field at several 
densities in the range 1.4 to 3xl0 10 cm~ 2 . The analy- 
sis uses the Dingle formula, AR/Rq ~ £/sinh£, where 
AR/Rq is the normalized amplitude of the SdH oscil- 
lations, £ = 2Tr 2 kBT/hu! c and uj c = eB/m*. Fits of 
the Dingle formula to the data are shown in Fig. 4, to- 
gether with the values of m* that provide the best least- 
squares fit. Within the experimental error (7%), m* is 
the same as the band effective mass of electrons in GaAs, 
mi, = 0.067m e , where m e is the free electron mass. For 
simplicity, in determining g* we have used the band ef- 
fective mass. 

In Fig. 5 we plot the values of g* , normalized to the 
GaAs band g-factor, as a function of density and r s [ p"5| . 
The open and closed symbols denote the two methods 
used to determine Bp\ their overlap confirms that the 
onset of exponential regime of the MR marks the full spin 
polarization. For the lowest density traces, we have data 
available only at 0.3K. However, T dependence of the 
in-plane MR, measured at n = 1.4 x 10 10 cm~ 2 , showed 
only a small variation of Bp with T in the range 25mK to 
0.5K. The larger error bars for the lowest density points 
in Fig. 5 take into account such variation with T. 

Figure 5 data highlight the main finding of our work: 
while the measured values of g* are three to six times 
larger than the GaAs band g-factor, this enhancement 
decreases with increasing r s . This trend is unexpected 
and at variance with theoretical calculations which pre- 
dict that, for an ideal 2DES, g* should increase as r s 
increases and the interaction becomes more important 
J|, ||. While we do not have an explanation for this 
discrepancy, we point out that at least three factors dis- 
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FIG. 4: Dingle plots of AR/R T vs. T. The density and 
the magnetic field at which the resistance was measured are 
indicated. The lines represent least-squares fits to the Dingle 
formula. For each dataset, m* obtained from the fit is shown. 



tinguish our samples from ideal 2D systems: finite layer 
thickness, disorder, and the spin-orbit interaction. 

For a 2DES with finite layer thickness, even in a single- 
electron picture, an in-plane magnetic field can affect the 
energy bands and therefore both the effective mass |lq| 
and the g-factor. Indeed, band calculations [O per- 
formed for the 2DES studied in our work reveal some 
enhancement of the effective mass and, to a lesser degree, 
a reduction of the g-factor with in-plane field. However, 
these modifications of the band parameters appear to be 
too small to explain the data of Fig. 5. 

The effect of disorder is more subtle. As the density is 
reduced, the disorder potential can play a more dominant 
role than the electron-electron interaction, and may lead 
to an inhomogeneous spatial distribution of the electrons 
in the sample. If the electrons become localized in the 
potential minima, it is possible that the behavior of the 
2DES reverts to that of a single-particle system and g* 
decreases. Such a scenario is consistent with the data of 
Fig. 5: for n = 0, g* /gb extrapolates to a value close 
to unity. On the other hand, we have made g* mea- 
surements on a sample with a mobility three times lower 
than the present specimen, and obtained similar results 
as those shown in Fig. 5. Moreover, our sample exhibits 
no signs of strong electron localization: it exhibits SdH 
oscillations, and its resistivity is less than h/e 2 . 

An examination of the results of several studies IO, 
|l2| , |l8| which have reported measurements of g* in 
Si-MOSFET 2D electrons provides further argument 
against disorder being responsible for the anomalous be- 
havior we observe in Fig. 5. These studies have gener- 
ally reported an enhancement of g* with increasing r s 
in much the same range of r s that we have examined. 
The mobility of our 2DES, however, is about a factor of 
10 larger than in Si-MOSFETs even though our densi- 
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FIG. 5: The effective g- factor normalized to the GaAs band 
g-factor (\gt\ — 0.44) vs. n or r s . The symbols denote the 
two different ways of measuring the full spin polarization field, 
Bp: from SdH oscillations in constant field (closed symbols) 
or from the onset of the exponential behavior of the magne- 
toresistance (open symbols). The dashed line is a linear fit to 
the lower density data points. On the right scale we indicate 
the corresponding values of AE/Ef, as defined in the text. 



ties are about 10 times smaller, implying lower disorder 
in our samples. If disorder were the main culprit, one 
would expect that Si-MOSFET data should also show a 
decreasing g* with increasing r s |19(| . 

Another important observation is worth describing. As 
we mentioned before, the in-plane MR of GaAs 2D holes 
was recently reported H, 0, p(l, and a relation between 
the MR behavior and full spin polarization, similar to 
that shown in Fig. 2, was established M. If we convert 



the magnetic fields above which the 2D holes are fully 
spin polarized to an effective g* , we find a trend very 
similar to the one seen in Fig. 5: g* is enhanced with 
respect to its band value but it decreases with increasing 
r s . We emphasize that the GaAs 2D holes too have very 
high mobilities and should contain low disorder. We con- 
clude that, while we cannot rule out the role of disorder, 
the unexpected behavior we observe in Fig. 5 appears to 
be intrinsic to low disorder GaAs 2D electrons and holes. 

The spin-orbit interaction, present in both GaAs 2D 
electron and hole systems, but nearly absent in Si- 
MOSFETs, may also play a role here. Since gb is sig- 
nificantly influenced by the spin-orbit interaction, one 
may expect that this interaction should modify g* in a 
many-body picture also. It is not clear, however, how the 
spin-orbit interaction would explain the trend in Fig. 5. 

Finally, in Fig. 5 we provide a measure of the inter- 
action energy for the spin polarization of our 2DES. We 
introduce the enhancement energy AE, defined as the 
difference between the band Fermi and Zeeman energies, 
both evaluated at the in-plane field where the 2DES be- 
comes fully spin polarized, i.e. AE = (2nh /mb)n — 
gb/J-BBp. The energy AE, measured in units of the Fermi 
energy is then simply equal to (1 — gb/g*). This quantity, 
which is indicated on the right scale of Fig. 5, reiterates 
the main result of our study. The measured AE/Ep de- 
creases with increasing r s , while theoretical calculations 
j|, || predict the opposite: AE/Ep should be zero at 
r s = 0, increase monotonically with r s , and reach unity 
for r s larger than a critical value of ~ 30. 
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